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The Boltzmann distribution (the most probable distribution) is one of the most important con¬ 
cepts used in physics, chemistry and biology. Suppose we put the system initially in one of the 
less probable state then the system will find the most probable state by a random search among all 
possible energy distributions and thus can take long time depending on the size of the system. In 
the following, simple analysis using our simple model shows that a small and physically reasonable 
energy bias against locally unfavorable energy distribution, of the order of a few kT, can reduce the 
time-scale of the process by a significant size. 


The Boltzmann distribution is one of the most important concepts used in physics, chemistry and biology [H-Q- A 
basic statement of the Boltzmann distribution in our context is, in a system of particles, the probability of a particle 
being in a state with energy e is given by 

P(e0 = (1) 

where g{ei) is the degeneracy of energy and /3 is a constant which determines the average energy. Also, /3“^ = kT, 
where k is Boltzmanns constant. In this paper, we assume the system is composed of a fixed number of particles 
denoted by ritot- Each particle is described by one or more one-particle quantum numbers, i, and a one-particle 
quantum energy, e^, a function of the quantum numbers. Particles need not be identical, but the total energy of the 
system is given, to a good approximation, by 

^total — ^ ^ niCi ( 2 ) 

i 

where rii is the number of particles with energy e^. We also assume the particles to be in dynamic equilibrium with 
sufficiently rapid exchange of energy through interactions between particles. For generality, we note that the quantum 
numbers, i, need not be one-particle quantum numbers; they could be collective excitations, e.g., phonons or normal 
modes. Suppose we put the system initially in one of the less probable state then the system will find the most 
probable state by a random search among all possible energy distributions and thus can take an enormously long 
time. In the following simple analysis using our simple model shows that a small and physically reasonable energy 
bias against locally unfavorable energy distribution, of the order of a few kT, can reduce the time-scale of the process 
by a significant size. We start with ntot number of particles and infinite number of energy levels, but only finite 
number of energy levels (say N) will be occupied. If rii is the number of particles having energy e^, then this number 
rii can be characterized as ’’correct” or ’’incorrect”, based on the final or perfect value Up (correct if rii = np, incorrect 
if rii ^ i^p). Correct number of particles are labeled as c, and incorrect number of particles are labeled i. A typical 
distribution of energy in terms number of particles in each of the total N energy levels can be expressed by a list of 

length N as ciciicicicciicc. It is obvious that the final or completely correct distribution is the one consisting of 

all c’s and no Fs. Now the problem is, starting from any arbitrary distribution of correct and incorrect number of 
particles at different energy levels, find how long it takes to get to the perfect distribution of energy for the hrst time 
? We assume that a correct value of number of particles at different energy level can become incorrect (c —>■ i) with 
the rate kd and an incorrect value of number of particles at different energy levels can become correct (z —>■ c) with 
the rate kic and that these changes occur entirely independently. As a result, the number (s) of incorrect number 
of particles at different energy levels changes in time. The first-passage time to the final state is the time required, 
starting from some arbitrary initial state s, to reach at s = 0 for the first time. The mean first-passage time r(s) 
is the average of this required time over all ways of getting from s to s = 0. We start with s incorrect number of 
particles in different energy levels, so the correct number of particles at different energy levels is n — s. The rate at 
which s —>■ s-|-1, is the correct number number of particles times the rate kd of changing a correct number of particles 
to an incorrect number of particles. 

Rate = kct{ri - s). ( 3 ) 

Similarly, the rate at which s —>■ s — I is the incorrect number of particles times the rate kic of changing an incorrect 
number of particles to a correct number of particles. 


(s—>8—1) 


Rate 


— kicS. 


( 4 ) 
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The probability that there are s incorrect number of particles at time t is denoted by P{s,t). This changes by gains 
from s — 1 and s + 1 and losses to s — 1 and s + 1. The master equation is 

= {n- s + l)kctP{s -l,t) + {s + l)fcjc(s + 1, t) - (n - s)kctP{s, t) - skicP{s, t), (5) 

with the assumption that P{—l,t) and P{N + l,t) are both equal to 0. Now we will use the standard procedure for 
using a master equation to find mean first-passage times. We start with the differential equations for P in matrix 
form as 


=^w{s,s')P{s',t). (6) 

s' 

Impose an absorbing boundary condition at S' = 0, so that only the states S = 1 to are involved. Then the 
fundamental equation that determines the mean first passage times is 

^t(so)w(so,s) =-1,; for all s. (7) 

So 


So that 


skic [r(s - 1 ) - t ( s )] + {n- s)kci [ t(s - 1 ) - t ( s )] = - 1 . ( 8 ) 

for all s between 1 and N. It is understandable that r(0) must vanish and t{N + 1) is never required. This determines 
all the other r(s). 


t{s) = -^{l + keq)^keq [ dy- - —{l + keqV) ^ \ (9) 

where keq = kd/kic- For large TV, the integral is dominated severely by the contribution from small values of y. It is 
very weakly dependent on s. Its asymptotic form for large N is given by 

= 1^(1 + [1 + ^'-(Nkeq)-^ + 2\iNkeq)-^ + .] ( 10 ) 

The s-dependent parts of r are in general negligibly small in comparison with the other terms, so the approximate 
result is given by 

t{s) Ki (1 + kci/kic)n. ( 11 ) 

flkci 

This asymptotic approximation is not valid for very small values of kd ■ The time t is now independent of the starting 
value of s, that is interesting because even if the starting configuration is close to final one, there is a significant 
probability that it will wander further away before reaching s = 0. The mean first-passage time for a fully biased 
search, where the change c —>■ f is not allowed so that kd =0, is given by 


r(s) 


he ^ j 


( 12 ) 


In this limit, t is independent of N and has a logarithmic dependence on s. Up to this point, the energy distribution 
was characterized only by N and the two rate constants (kd and kic). Now, it will be interesting if we understand 
the deeper meaning of the ratio kd/kic- The kinetics for number of particles in any energy level is given by 

with [c] -I- [z] = 1. (13) 


i^eq 
kic 


[c] — kd [c] -f kic [f ] 


So the equilibrium constant is defined as 


(14) 
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So that 


leg 


1 + keg 


and 


lleg — 


1 + keq 


(15) 


It is easy to understand the concept of equilibrium constant keg from statistical physics. Suppose the free energy for 
correct number of particles in the i-th energy level is Gc, and free energy for incorrect number of particles in the 
same energy level is given by Gi = Gc + U. So U is an free energy increament for having an incorrect number of 
particles, U is definitely a positive number but it magnitude depends on the extent of ’incorrectness’. Then by using 
the equilibrium statistical thermodynamics, one gets 


. _ ^ _ -U/kT 

l^eg — , — e 


(16) 


So that kci/kic is small, if U is larger than kT and r can become much smaller. Also at higher temperature the system 
will attain Boltzmann distribution faster. 
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